Coarsening Dynamics of a Quasi One-dimensional Driven Lattice Gas 
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We study domain growth properties of two species of particles executing biased diffusion on a 
half-filled square lattice, consisting of just two lanes. Driven in opposite directions by an external 
"electric" field, the particles form clusters due to steric hindrance. While strictly one-dimensional 
systems remain disordered, clusters in our "quasi ID" case grow until only a single macroscopic 
cluster survives. In the coarsening regime, the average cluster size increases ~ t 6 , significantly 
faster than in purely diffusion-controlled systems. Remarkably, however, the cluster size distribution 
displays dynamic scaling, following a form consistent with a diffusion-limited growth mechanism. 
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Introduction. The dynamics of a system undergoing 
phase segregation when quenched below the transition 
temperature has been of interest for many years. Start- 
ing from a homogeneous state, domains of the co-existing 
phases form and grow. Nearly all existing studies of such 
coarsening processes are devoted to systems subjected to 
dynamical rules which eventually take them to equilib- 
rium states Our interest is coarsening in systems 
which are evolving toward nonequilibrium steady states 
(NESS): The underlying dynamics of our systems vio- 
lates detailed balance and time reversal invariance. In 
particular, our main focus will be the remarkable char- 
acteristic of dynamic scaling. In this letter, we report 
the presence of this behavior in a simple model of biased 
diffusion of two species. In stark contrast to similar phe- 
nomena in other systems where the growth exponent is 
typically 1/2 or 1/3 [§-[§, we find it to be at least 0.6! 
This large, anomalous value rules out diffusion-controlled 
growth mechanisms. Surprisingly, however, the scaling 
function for the cluster-size distribution resembles that 
of diffusion-dominated growth ||,^| quite well. Since our 
system is essentially one-dimensional, we may also ex- 
ploit a reduced description, in terms of coalescing ran- 
dom walkers on a line. To account for the faster growth, 
we introduce interactions between neighboring walkers. 
Apart from small deviations, the "cluster-size" distribu- 
tions in this picture are in good agreement with the data. 
In the following, we briefly describe the model, present 
the data from Monte Carlo (MC) simulations, discuss the 
coalescence picture, and conclude with some open ques- 
tions. More details will be left to ||. 

Biased Diffusion of Two Species. Motivated by the 
physics of fast ionic conductors Q, Katz, et.al. || gener- 
alized the well-known Ising lattice gas to include a uni- 
form bias affecting all particles. With periodic boundary 
conditions, the system settles into a NESS rather than a 
simple equilibrium state. In addition to a non-trivial par- 
ticle current, many novel collective phenomena emerge, 
some of which remain unexplained |9|]. A natural ex- 
tension of such "driven diffusive systems" involves two 
species. When these are driven in opposite directions, 
the system displays a phase transition in d — 2, even if 



the only interparticle interaction is an excluded volume 
constraint [|l0|,^TJ . In d — 1 , however, one can prove fj^| 
that phase transitions cannot occur. Since such a system 
consists of particles on a line, it may be used to model 
traffic flow p3[ , and we will use the term "single lane" 
for this d = 1 case. Remarkably, when a second "lane" is 
introduced, the system again displays a phase transition 
|Q ! We briefly describe this study. 

A fully periodic 2 x L square lattice is randomly filled 
with equal numbers of two types of particles (labeled by 
"charge" + or — ), subject to an excluded volume con- 
straint. When not driven, the particles diffuse accord- 
ing to the following rules. Two nearest neighbor (NN) 
sites are chosen at random, and attempts to exchange 
their contents are (a) always allowed for particle-hole 
pairs and (b) accepted with probability 7 for particle- 
particle pairs. Under this dynamics, the two species are 
just diffusing randomly. The system obviously remains 
homogeneous, with trivial collective properties. Next, we 
impose an external "field" which drives the two species 
aroung the "ring" in opposite directions, by simply for- 
bidding all exchanges which result in, say, +/— particles 
moving clockwise/anticlockwise. In Q|, 7 is fixed at 0.1 
while systems up to L — 10 are evolved for as many as 
4 x 10 6 MC steps (1 MCS corresponds to 2L exchange 
attempts). For early times (10-20 MCS), small block- 
ages form everywhere, due to the mutual obstruction of 
the opposing species. These clusters coarsen until, at 
late stages when NESS is essentially established, only 
a single macroscopic cluster remains. The characteris- 
tics of this cluster are notable, (a) It contains almost 
no holes, (b) It exhibits a non-trivial "charge" profile, 
(c) Its size (number of particles connected as NN's) is 
gaussian distributed around 0.94L, so that its "length" is 
about L/2. (d) A low-density region spans the remainder 
of the system, consisting of particles ( "travellers" ) which 
leak out of the cluster at one end and later rejoin it at 
the other. By contrast, for the "single lane" case, the 
system remains homogeneous, and the cluster size distri- 
bution is known exactly fl^ | to decay exponentially (with 
weak 1/L corrections). In this sense, we believe that it 
is justifiable to use the terms "long range order" and 
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"phase segregation" to characterize the two-lane system 
in steady state. Of course, if the overall particle density 
is reduced sufficiently, there must be a phase transition to 
re-establish homogeneity. In an effort to understand why 
a two-lane system orders while the single-lane version re- 
mains disordered, we turn to a study of how long-range 
order builds up in the latter case. 

Simulations, Characterizations, and Results. To study 
the time dependence of our stochastic system, we carry 
out many runs with the same random initial distribution 
of particles. Typically, data from 100 — 1000 runs are 
used to build histograms. In addition, with finite size 
effects in mind, we exploit a range of system sizes, from 
2 x 100 to 2 x 1000. At half filling, local blockages ap- 
pear everywhere soon after the run starts, a phenomenon 
reminiscent of homogeneous nucleation. Though small, 
these blockages share some characteristics of the termi- 
nal cluster, having few holes inside and growing through 
accretion of particles from the inter-cluster region. In a 
relatively short time (~100 MCS), a rough balance be- 
tween accretion and loss is established so that both the 
density of travellers and their total number remain rel- 
atively constant for the rest of the run. The associated 
values, ~ 5% and ~ 0.06L respectively, can be under- 
stood qualitatively Thus, we may characterize this 
initial stage as one in which local densities quickly reach 
their (approximate) final values. For later, intermediate 
times (the "growth" regime), the larger clusters grow at 
the expense of smaller ones, a mechanism reminiscent of 
Lifshitz-Slyosov growth. However, the similarity is de- 
ceptive, as we will show next. 

To characterize domain growth in our model, we mon- 
itor the cluster size distribution since it carries very de- 
tailed information. The most natural definition of the 
size of a cluster, s, is just the number of particles (regard- 
less of their charge) connected via NN bonds. We collect 
data at various times, t, to construct histograms of such 
cluster sizes: p(s, t). Since this distribution is sensitive to 
both the growing clusters and the travellers, we refer to it 
as the microscopic cluster distribution. Clearly, p is not 
conserved, decreasing when clusters merge, etc. Instead, 
we consider p(s,t) = sp(s,t), which, counting just the 
total number of particles, is conserved. Known as the 
"residence distribution" in percolation studies, p(s,t) is 
proportional to the probability of any site belonging to a 
cluster of size s at time t. For the 2 x L case, p displays 
two peaks from an early time, one at s = 1 (followed by a 
sharp drop to nearly zero) and another moving to higher 
s with t. As mentioned above, the first is relatively con- 
stant, so that we focus only on the behavior of the second 
peak. To this end, we devise a coarse-grained (CG) de- 
scription of the clusters, which has the added advantage 
of tolerating the occasional hole in a cluster. For any 
configuration on the 2 x L lattice, we construct another 
one, with occupation numbers or Ion a line of L sites, 
as follows. At each site i, we assign if there are five 



or less particles in the ten sites around the i th column of 
the original lattice; and 1 otherwise (cf. Fig. 1.) 
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FIG. 1. Illustration of the coarse-graining procedure for a 
configuration in a L = 12 periodic lattice. An example of 
the microscopic configuration is shown with ± particles and 
holes. The CG version shows only {0, l}'s. 

In the CG configuration, a cluster is simply a consecutive 
sequence of l's and its "size" (£) is just the length of this 
string. Histograms for these cluster sizes, p(£,t), can be 
compiled easily. Note that clusters in the traveller region 
are usually quite small (s < 5), so that p(£,t) has only 
one peak, associated with the growing clusters. Again, 
we consider the residence distribution: p(£,t) = £p(£,t), 
from which a natural definition of an average cluster 

length arises: £ c (f) =T,^-p(^)/HtP{^t)- In Fi S- 2 > 
we present a log-log plot of this quantity for various L 
(solid points) . Clearly, there is good data collapse in the 
growth regimes. For the two smaller systems, our runs 
are long enough for saturation {\.m\ t ^,oo £c{t) — 0.47L) 
to occur. More significantly, the effective power in the 
growth (i.e., din tc /dint) appears to be still increas- 
ing for the largest system at the latest times. In the 
figure, we provide two straight lines, corresponding to 
t x l 2 and t 2 / 3 , to guide the eye. From these lines, we 
estimate that the effective exponent is about 0.6 at the 
latest stage. Unless a mysterious slow-down sets in at 
even later times, we are confident that the power of 1/2 
is ruled out. Finally, we turn to the distributions and, 
to check for dynamic scaling, construct the normalized 
scaled form, / (x, t) = £ c (*) p(£, t) with x=£/£ c (t). As 
we see in Fig. 3a, / (x, t) is largely independent of t (and 
L) in the growth regime. Thus, we conclude that, de- 
spite the possible absence of pure power law growth, the 
distribution displays, remarkably, dynamic scaling. Fur- 
ther, the scaled distribution is well approximated by the 
solid line in Fig. 3a, a distribution (cf. Eqn. (1) below) 
well known from diffusion- limited coalescence [ 3 Q . This 
suggests casting our MC dynamics in terms of an effective 
model, involving coalescing random walkers (RW) . 

Coalescing random walkers in d = 1. Since the density 
of travellers is relatively constant, a good macroscopic 
picture focuses on the sizes (lengths) of the clusters only, 
with the main dynamics being particle exchange between 
neighboring clusters. To describe this sequence of clus- 
ters, points ("walkers") are placed on a ring of N sites, 
with spacings corresponding to the cluster-lengths [pi. 
Particle exchanges translate into walkers taking steps. 
When two walkers meet, they "coalesce," representing 
the disappearance of a cluster. The number of walkers, 
n, maps into the number of clusters, while N represents 
the sum of all cluster lengths. At the end, there is only 
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one walker, i.e., one macroscopic cluster. For simplicity, where C is an amplitude giving the best fit to the data. 



we set N = L/2 which well approximates the observed 
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FIG. 2. Coarse-grained MC (solid symbols) and interacting 
random walk (open symbols) data for the average cluster size 
vs time, for L = 10000, 1000 and 500. 

If the random walkers were free (apart from coales- 
cence), the problem is exactly soluble in d = 1 pjl. De- 
noted by Po(£,t), the (unnormalized) frequency of find- 
ing two adjacent walkers separated by a distance £, at 
time t, satisfies a diffusion equation in the continuum 
limit. With absorbing boundary condition at £ = to 
model coalescence, the solution is standard. To compare 
with the MC data, consider the normalized "residence" 
distribution, Pq(£, t) = £Po(£,t). Its first moment, £ , 
should correspond to tc- Clearly, £q grows diffusively: 
£o(t) oc t 1 / 2 . Further, dynamic scaling prevails (in the 
limit t — > oo, £ — > oo at fixed x = £/£o(t)) : 

32 4 
£ (t)P (£,t) = F (x) = -^z 2 exp( x 2 ) , (1) 

which serves as a good benchmark (solid line in Figs. 3a 
and b) for scaling plots of the data. Unfortunately, the 
growth exponent (1/2) of this simple model is unequiv- 
ocally ruled out by our data (Fig. 2). Its deficiencies 
can be traced, by monitoring the evolution of our clus- 
ters, to the rapid, systematic disintegration of small ones. 
To model this bias towards larger clusters, we introduce 
a model of interacting random walkers (IRWs), where 
neighboring walkers experience stronger attraction with 
decreasing separation. Specifically, consider a walker and 
its two nearest neighbors. Letting £r (£l) be the distance 
to its right (left) neighbor, and q R , q L , q s (= 1-qR- qh) 
be the probabilities for moving to the right, left, and 
staying, respectively, we choose: 
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FIG. 3. Scaled residence distributions for CG MC, f(x), 
and IRW, F (x). For comparison, the "benchmark," Fo, is 
shown as a solid line in both plots. In (a), ■ A and • 
are data from a 2 x 2000 lattice at t = 1096, 2981, 8103, and 
2207, respectively. Open symbols are data from 2 x 500 at 
t = 403, 1096, 2981, and 8103. In (b), this set of solid (open) 
symbols correspond to TV = 1000 (500), at scaled times. 

Some comments are in order, (i) To motivate the in- 
dependence, we invoke the mean-field theory for the two- 
species model pd| ] . Considering the continuity equations 
for mass and charge density for the traveller region and 
the interior of a cluster separately, and matching the two 
regions via a shock condition we find the steady- 
state charge current through a cluster of length £ » 1 to 
be j{£) = (7/2) [1 + 0(£~ 2 )] §]. Symmetry dictates that 
the mass currents of the +/— species are equal and oppo- 
site. Thus, for neighboring clusters with in and £l, say, 
the imbalance in the currents would be ~ {£/^ — ^ 2 ) . 
In turn, this current difference leads to one cluster gain- 
ing at the expense of its neighbor, (ii) To account only 
for a bias, qs = would suffice. Our choice, allowing a 
walker between small clusters to be more mobile, mod- 
els the faster dynamics observed in that case. (Hi) If C 
were a simple constant, the behavior at late times (with 
fewer but larger clusters) would again be dominated by 
^-independent rates, leading to t 1 !" 1 growth. As a rem- 
edy, we must supply a scale, namely C, against which 
to measure "small" clusters. Thus, C is chosen to de- 
pend on the average size of clusters, i.e., N/n. Finding 
C ~ (0.06N/n) to fit best, we provide the following 
heuristic argument in its favor. On the average, an £- 
cluster will find itself sandwiched between two sections of 
travellers with length N/n. When £ <C N/n, the overall 
density in that region will be too low to support inhomo- 
geneities and the cluster should decay quickly. As for the 
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value 0.06, we note that this is also the traveller density, 
but can offer no insight into this coincidence. 

This choice for C renders our rates time-dependent 
(through nit)), so that the exact solution of this dynam- 
ics remains elusive. Fortunately, our IRW is easily sim- 
ulated. At time t = 0, walkers are distributed randomly 
over the ring, with an initial density of 0.4, and then 
updated according to (|J). We measure the normalized 
residence distribution P(£,t) and its first moment, £i(t). 
Adjusting a single scale parameter to match t to t, we 
plot (in Fig. 2) both lc and li, for a range of system 
sizes. The success of our effective model is quite stag- 
gering: the two measures for £ trace out the same curve 
for over four decades. The small differences observed in 
the steady state values are easily understood: while £i(t) 
approaches a final value of L/2, the corresponding limit 
of ic{t) is only ~ 0.47L. Finding P{£,t) to be qualita- 
tively similar to its counterpart from MC, we construct 
the scaling plot: F(x) = h(t)P{l,t) vsie £/£i(t) (Fig. 
3b). Within statistical noise, the data collapse is quite 
satisfactory. Most remarkably, this F(x) also closely fol- 
lows F (x) of the free RW! 

There are, however, small deviations from scaling, re- 
vealed only by close inspection. We first note that, in 
both cases, the late-time data tend to deviate systemat- 
ically from Fq(x). In Fig. 3a, the data lie slightly above 
Fq for x < 0.5, and slightly below Fq for x near its maxi- 
mum. In Fig. 3b, the trend is reversed. For x > 1, there 
are no perceptible deviations in either case. While we 
cannot offer any rigorous understanding, it is conceivable 
that the lengthening domains of travellers - which are ne- 
glected in the IRW dynamics - play a role here. Another 
possible mechanism is described in [ |l6| . Much more re- 
markable, however, is how close both / and F follow Eqn. 
(|l|), i.e., the free RW scaling function, even though nei- 
ther £c and £i are consistent with the diffusion equation 
leading to Fq ! It appears as if the complexity of our full 
dynamics can be absorbed into a single "renormalization" 
of time, t (x T a , such that simple diffusion-controlled co- 
alescence re-emerges in r. Matching t x / z ~ £ ~ t 1 / 2 then 
requires a = z/2. Clearly, formalizing this scenario poses 
a serious theoretical challenge. 

Summary and Outlook. To conclude, we summarize 
our findings and end with some speculations. We con- 
sider a lattice gas with two species of particles driven in 
opposite directions on a "single-lane road" vs. a "two- 
lane highway." Allowing a small amount of particle- 
particle exchange, both systems are ergodic but display 
drastically different final states: While the one-lane sys- 
tem remains homogeneous, the two-lane case exhibits a 
macroscopic cluster which scales with system size. In an 
effort to probe the surprisingly different outcomes, we 
study the growth of clusters in the two-lane case in some 
detail, measuring the average cluster size, £(t), and the 
full residence distribution. We discover that £{t) grows 
much faster than typical domains in similar models. At 



late times, the growth law may be characterized by an 
effective exponent of about 0.6. We end with a specu- 
lative note. In ordinary Lifshitz Slyosov growth, which 
leads to the power 1/3, the mechanism for transport (of 
material from smaller clusters to larger ones) is diffusion. 
Here, particles suffer biased diffusion, so that the motion 
is ballistic instead. Can the difference of a factor of two 
between these different mechanisms transform the power 
1/3 to a 2 /3? Needless to say, to arrive at definitive con- 
clusions, we need not only a deeper analytic understand- 
ing but also a better Monte Carlo study of this simple, 
yet intriguing, model system. 
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